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WORLDWIDE AFFAIRS 


BRIEFS 


CHINESE, IRAQI COOPERATION ACCORD--Baghdad, October 10 (XINHUA)--An agree- 
ment on information cooperation between the new China (XINHUA) news agency 
and Iraqi News Agency was signed here this morning. It was signed by 

Hou Yefeng, Chinese ambassador to Iraq, on behalf of the XINHUA News Agency, 
and Muadh Abdal-Rahim, acting director general of the Iraqi News Agency. 
[Text] [Beijing XINHUA in English 1200 GMT 10 Oct 79 OW] 


PRC DELEGATION TO BELGIUM--Paris, October 6 (XINHUA)--The Chinese Posts 
and Telecommunications delegation led by Wang Zigang, minister of posts 


and telecommunications, left here for Brussels today after a friendly 
visit to Franee, On October 4, a protocol of the agreement on cooperation 
between the two countries in the field of telecommunications was sipned 


by Minister Wang Zigeng and secretary of state for Posts and Telecommunica- 
tions of France Norbert Segard. [Excerpt] [Beijing XINHUA in English 
0248 CMT 7 Oct 79 OW] 


"INA,’ "NCNA* SIGN AGREEMENT--Baghdad, 10 October--INA and the Chinese News 
Agency, NCNA, signed here today an agreement on cooperation and the ex- 
change of news and correspondents. The agreement stipulates the exchange 
of news bulletins and cooperation in rendering free services and assistance 
to the correspondents of the two news agencies in the two countries. 

Ma'adh ‘Abd al-Rahim, acting chairman of INA board of directors, signed 

the apreement for TNA, and PRC Ambassador to Iraq Yefeng Hou signed for 
NCNA. [Text] [Baghdad INA in Arabic 1030 GMT 10 Oct 79 JN] 


LIBYAN RADIO STATION--Tripoli, Oct 4, JAMAHIRIYAH News Agency--Mr Muhammad 
al-Shuwaydi, the JAMAHIRIYAH secretary of the General People's Committee 

for Information, signed a contract with an Italian company to install a 
radio station in Tobruk. The information secretary had imported the 

station equipment and instruments. This radio station will transmit to 

many parts of the Arab world. The project will cost 1,648,356 Libyan dinars 


and will take one year for completion. [Text] [Tripoli JANA in English 
1913 GMT 4 Oct 79 LD] 





UCANDA-USSR INFORMATION AGREEMENT--The minister of information and broad- 
casting, Dr Albert Picho Owiny, has met the ambassador of the Soviet Union, 
Mr Sergey Bukin, with whom he discussed the possibility of renewing the 
news and programs exchange agreement due to expire this vear, The news 
agreement was signed in 1977 between the Ministry of Information and 
Broadcasting and the Soviet Union News Agency. TASS. The agreement which 
is renewable after every 2 years allows the Ministry of Information to 
receive and use news and information transmitted by TASS in English and 
the TASS correspondent in Uganda to file and to receive news from the 
Ministry of Information and Broadcasting in Uganda, Dr Owiny and Mr Bukin 
agreed in principle to renew the agreement, The minister and the 
ambassador also agreed to meet again for further discussions related to 
the agreement. The meeting was attended by the permanent secretary, 
Ministry of Information and Broadcasting, Mr Nyakimwe; the director of 
information, Mr Nathan (Akenu) and two officials from the Soviet Embassv. 
|Yext| [Kampala Domestic Service in English 1700 GMT 3 Oct 79 LD/EA] 
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JAPAN 


BRIEFS 


INCREASED TELEPHONE USE--Geneva 20 Sept KYODO--Telephones will be 
installed for one of every two persons in Japan by the end of the 
century, the chairman of the Nippon Electric Co predicted Wednesday. 
Speaking to the Third World Telecommunication Forum, Dr Koji Kobayashi 
said Japan operated 40.8 telephones for every 100 persons in 1977. The 
United States operated 71.8 and Switzerland operated 63.8 telephone sets 
per 100 persons, he said. Kobayashi said the 20 percent increase in 
telephone use will be aided by the large integration of solid state 
equipment, practical use of opto-electronics and the combined use of 
computers and communication equipment. He also said more diversified 
use of satellites will aid telephone communication. [Excerpt] 

[Tokyo KYODO in English no time given 20 Sep 79 OW] 
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MALAYSIA 


BRIEFS 


SATELLITE RECEIVING STATION--Petaling Jaya, Peninsular Malavsia, 8 Oct 
(BERNAMA)--A $1.5 (US$600,000) million satellite receiving station will 

be set up in Kuala Lumpur next year, Deputy Transport Minister Dr Goh 
Cheng Teik said on 8 October. He said the existing tracking station in 
Kuantan was no longer adequate. With the installation of the new station, 
the meteorological department should be able to receive the clear pictures 
of the atmosphere taken by Japanese and American satellites orbiting the 
carth. The informatio, obtained would help to identify tropical dis- 
turbances in the South China Sea and elsewhere before they hit the Molavsian 
coastlines. “We will not only be able to provide early flood warnings 

to the population but also advise planes, ships and patrol craft on sea 
and weather conditions," Dr Goh said. "We will also be able to alert the 
oil rigs operating off Trengganu, Sabah and Sarawak." [Excerpt] [Kuala 
Lumpur BERNAMA in English 0816 GMT 8 Oct 79 BK] 
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THATLAND 


TELEPHONE CHIEF DISCUSSES COMPUTERIZED TELEPHONE SYSTEM, OTHER POLICIES 
Bangkok BAN MUANG in Thai 22 Jul 79 pp 1, 16 
/Arcicle: “Telephone Improvement s"/ 


/Text/ Mr_Boonchu Phianphanit, director of the Telephone Organization of 
Thailand /TOT/ talked to BAN MUANG reporters at the Ministry of Communica- 
tions and answered questions regarding substantial increases in the price 
of oil which has forced the majority of businessmen to carry out their 
business transactions by telephone. When asked what service the TOT was 
prepared to provide the people in alieviating the problem, Mr Boonchu 

said that currently 1.3 million calls are placed each day from 300,000 
telephone numbers in Bangkok and the vincinity and that the businessmen 
can make long distance calls to carry owt their business transactions, but 
TOT cannot help them with transportation problems. 


A BAN MUANG reporter mentioned that operators need to improve their 
efficiency because they often connect calls with wrong numbers even though 
callers dial the right ones. Mr Boonchu said the problem has not been 
resolved, and that it is a matter of personnel efficiency. We said that 
once caught misconnecting the numbers the operator's salary would be 
docked, and in the past there had been several firings; such punishments 
should therefore help the matter. However, if one dials the right number 
and gets another the customer must inform chief of the relay center in the 
area se that the telephone circuits can be corrected at the central relay 
center. 





As to the question cf whether it is true that members of the House of 
kepresentative request several telephones but later lease them to merchants, 
Mr Boonchu said that each member of the House must submit a written request 
approved by the Secretary General of Parliament to TOT to have one phone 
installed in his provincial residence and another at this Bangkok residence 
only, and about 10 percent of these representatives have submitted the 
requests. Regarding the allegation that the representatives lease their 
phones for profit, Mr Boonchu siad, “They are honorable members. If they 
request telephones they shov'd use them properly since TOT has no way of 
checking on them. We provide phones to whoever carries the letter from 

the Secretary General of the Parliament. The media has to work as our 

eyes and ears if there is any leasing going on." 


5 











The BAN MUANC reporter further asked when the telephone system will be 
computerized, Mr Boonchu responded that within a month or so TOT would bid 
for 100,000 numbers ~it*‘a a budget of 10 billion baht. TOT has received 

a loan of 1,800 millic« irom the World Bank, 1,200 million from the Bank of 
Japan, 3,600 million from varicus commercial banks, the balance is in the 
form of loans from the sale of bonds. The actual costs of installation are 
not yet known until the bidding is accomplished, they may be bighez or 

the same as the previous rates. As for automated long distance service, 

it will be initially operated between the provinces of Chiang Rai, Chiang 
Mai, Lampang, Phrae, Denchai /District] and Nan Province, and during March 
and April of next year to Bangkok, and the service is expected to be 
operated nationwide by the yea. < “80. 


Mr Boonchu further disclosed thet ‘fore the end of the year, \»,%, public 
telephones will be operated at 4 { sio.;2 Limit, and if a caller desires 
to continue, he must deposit anothe; coin, Many complaints have been 
received that some people would use the public phone for nearly half an 
hour at a time. However, home telephones present no problems so far. | 
TOT is trying to figure out the average use time and it appears that 
customers are using telephones at length, they must also have imposed a 
time limit. As for income, last year (1978) TOT earned 1,238.6 million 
which was 123.6 million more than 1977, and it certainly hopes to see an 
increase in 1980, Mr Boonchu said. 


9013 
cso: 5500 








LANDSAT PHOTOS--Mr Suwit Wiboonset, director of the National Resource 

Surveys Via Satellite Division disclosed that around the year 1981 

Thailand will have a station to receive satellite photos or photo signals 

directly fro the United States’ LANDSAT which provides data on natural | 
resource surveys. Mr Suwit also said that this would be economically 
beneficial since the data can be received in a timely fashion. Now it 

sometimes takes 6 months for Thailand to receive photos from NASA. The 

future station, according to Mr Suwit, is capable of relaying natural 

resource photos, as well as changes in the land every 18 days. Mr Swit 

further disclosed that besides Thailand, other Asian nations will also 

benefit from this station since LANDSAT is capable of taking pictures over 

a wide area, from China to Bangladesh, Sri Lanka, Burma, Thailand and 

other ASEAN countries, and these countries can contact Thailand if they 

are interested in_sharing the photos. /Text/ /Bangkok SIAM RAT in Thai 

23 Jul 79 p 1, 12/ 9013 


LONG DISTANCE CAPABILITY--The Telephone Organization of Thailand will 

expend its long distance service throughout the country without having to 
commect calls to the Switching boards by the year 1980. Mr Chalor Mancerin, 
Operating Director of TOT told SIAM RAT BUSINESS WEEKLY, regarding the 
expansion of the national long distance service that by the beginning of 

1980 the organization would open service to six northern provinces without 
having to pass throwgh the switchboard. The next project includes expanding 
the service to enalle customers to dial direct throughout the country, 

that is, callers can place a long distance call by merely dialing the 

desired number. “Mr Chalor further disclosed that bv the middle of next 
year, 1°80, customers in provincial areas should be able to call Bangkot 
numbers directly, and that by the end of 1980 telephone communication between 
hangkok and upcountry should be fully in operation. At present, customers 

in a certain province, such as Pattaya in Chonburi Province are already 
miking direct calls within their province. /Text/ /Bangkok SIAM RAT BUSINESS 
WEEKLY in Thai 30 Jul - 5 Aug 79 pp 1, 4/ 9013 
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CUBA 


BRIEFS 


SANTIAGO TELEVISION TRANSMITTER--A television signal retransmitting station 
taking black-and-white and color programs to the southwestern coast will be 
built in Santiago de Cuba Province, Guama Municipality. (Mariano Romano), 
director of the communications enterprise in Santiago de Cuba, reported 
that the civil construction will start soon in the place known as “Punta 
Tabacal" located some 4 km from the village of Chivirico. About one year 
ago experimental equipment was placed in this area to observe the behavior 
of television signals and the satisfactory results led to set up the in- 
stailation. [Text] [Havana Domestic Service in Spanish 2300 GMT 2 Oct 79 
FI] 
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LIBYA 


BRIEFS 


SEA CABLE, GROUND SATELLITE STATION--Tripoli--Last night, the secretary 
of the General People's Committee for Communications and Maritime Trans- 
port opened the sea cable system between Marseilles and Tripoli, and 
Benghazi and Tripoli, and the ground satellite station. The opening 
ceremony was also attended by the ambassadors of Japan and France and 

by the directors of French cable and wireless and the Japanese Nippon 
clectricity firm. The sea cable system between Tripoli and Benghazi 
will provide 900 multichannels for telephone and cabie communications 
between the two towns in addition to the simultaneous relay, both ways, 
of radio and colored television programs side-by-side the microwave net- 
work all along the Libyan coast. The sea cable between Tripoli and 
Marseilles, will provide telephone and cable communications between 

the Jamahiriyah and France and from there to the rest of the world. The 
sround satellite will provide telephone and cable communications and 
televisiopv programs between the Jamahiriyah and the states subscribing 
in the satellite communications. [Tripoli JANA in Arabic 0900 GMT 

20 Sep 79 LD] 
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FRANCE 


BRIEFS 


SIEMENS, THOMSON TELEPHONE CONTRACT--The German Siemens group and the French 
Thomson group have recently signed an important agreement with the Egyptian 
authorities. This agreement, amounting t> a total of 1.8 billion dollars 
(7.5 billion francs) deals with the renewal of the Egyptian telephone net- 
work and the installation of 500.000 new lines originating from "temporal" 
centers which will permit to triple the capacity of the network. The agree- 
ment is encompassed within the framework of a plan of 5 years" duration due 
to start in 1980. Siemens’ share in this contract will be about 1.1 billion 
dollars, with a large part of the work due to be financed by the West German 
aid development fund. Thomson's share will be about 450 million dollars. 
American firms will also take part in the operation with a share of 250 mil- 
lion dollars in connection with credits granted at favorable rates by 
Washington. [Text] [Paris LE MONDE in French 19 Sep 79 p 42] 8696 


TELECOMMUNICATIONS DIRECTORATE DECENTRALIZATION--The directorate of tele- 
communications of the national network (DTRN), presently located in Paris, 
is in the process of decentralization in the air-port zone of Blagnac (Hte 
Garonne) in 1981. Five hundred telecommunication employees are affected 
by this move. The directorate is in charge of the programming, studies, 
construction, utilization and maintenance of the interurban network for 
medium and long distances. Side by side with this Blagnac installation, 

5 operational directorates will be created in Lyon, Metz, Nantes, Paris 
and Toulouse. Agreed upon on 1 March 1977 by an interministerial committee 
of land management, the decentralization of DTRN had at first provoked a 
strike movement among its personnel who feared "exile in the province” 

(LE MONDE 2 March 1977). The lack of repeated trouble seems to demonstrate 
that consultations with trade union organizations allowed to lift prejudice 
about a decision necessary for land management. [Text] [Paris LE MONDE in 
French 19 Sep 79 p 31] 8696 
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NUMERICAL CALCULATION OF RAY OPTICAL PATHS, RELATIVE DELAYS IN 
GRADED-INDEX OPTICAL FIBERS 





Milan ALTA FREQUENZA in English Jul 79 pp 240-2458 


[Article by G. Cancellieri, P. Fantini, Fondazione G. Marconi, Pontecchio 
Marconi/ 


[fext/ Abstract . Using an intrinsic tern of variables, it is pos; ible to give a 
simple and evident description of ray optical paths for whichever circularly 
symmetric index profile, and to calculate ray relative delays, even for near- 
-optimum index profiles.,These numerical calculations can be very useful in 
order to explain the experimental results of 4 differential mode measurement, 
and to find some new applications of such 6 measurement. 





1. INTRODUCTION 


The variables which are usually assumed for cescri- 
bing ray opticel paths in a graded-index optical 
fiber are z, r. ¥ where z is the axial coordinate, 
r the radial coordinate, and » the azimuthal angle. 
Each ray optical path is function of some initial 
conditions which involve the direction of the ray. 
In this paper, we develop @ model of the rey prope- 
gation on the basis of an intrinsic tern of verie- 
bles r, @. 6. where ® is the angle between the lo- 
cal direction of the rey and the z-axis, and ¢ is 
the angle between the radius r and the projection 
of the ray on the trensverse plane. 

Such e description offers two main advantages: 
firstly, the initial condition can be fully expres- 
sed as (roe 4," 0.) without introducing any other 


variebles secondly, ali these variables ere perio- 
dic z-functions, therefore the description of one 
period is sufficient for a complete characteriza- 
tion of the ray optical path. As e direct conse- 
quence of the first property. this description seems 
to be particularly useful for differential mode me- 
asurement, in which the initial condition (rr. @., 
e) can be considered wither a launch condition 


[1,2]. or 6 detection condition [3,4]. 
Once the generic ray optical path over «6 period Zz 
is obtained, it is possible to evaluate the ray re- 


ll 




















lative Geley. ‘his new method for @ theoretical in- 
vestigstion of modal dispersion in ea graded-index 
fiber results as eccurste es the previously propo- 
sec ones. Its main feature is in the possibility of 
relating each computec relative delay to a well de- 
fined geometrical trajectory. From the theoretical 
point of view, this can simplify the study of the 
behaviour of fibers having whichever index profile. 
Moreover, this model can be considered as a ster- 
ting point for extending the investigation to-wa- 
velength Gependent index distributions; e future 
paper will be devoted to this eatension. From the 
experimental point of view, two apolications are 
actually under study. in the field of differential 
mode delay measurement. 

The first one starts from 6 very eccurete measu- 
rements of the index profile nir),. and consists of 
explaining the differences between theoretical and 
experimental relative delays. in terms of mode cou- 
pling and differential mode attenuetion {s}. The 
second one, starting from the knowledge of the in- 
dex profile only as a first approximation, consists 
of taking advantage from this model, in order to 
enhance the precision of the profile measurement. 
Using only meridional rays. this technique hes bees 
already proposed [2]. Taking into account also skew 
rays. to get a higher precision soems possible. 

We report here some exomplifying results. espe- 
cially for whet concerns near-optimum profiles, for 
which a very high precision is required. In sect. 2, 
the outline of this model wili be presentcc, devo- 

ing perticuler attention to helical rays. In Bect. 
3, the optical paths of some characteristic rays 
will be investigeted, to the point where it is pos- 
sible to obtain comparative plots of r, @. 6. as 
functions of z, over @ period Zz. Finally. sect. 4 
will be devoted to the relative delay calculaetion. 


2. PROPAGATION DESCRIPTION BY THE INTRINSIC TERN 
Tv. o. 5. 


For simplicity. we take into account 4 medium he- 

ving = <0, for O<r¢a, and an * 0, for roa. Letting 
n{o) = n. and nia) — ne end considering en initiel 
condition (rie ¢,” 0.) the ligth acceptance of the 


fiber is defined by the relation 


12 


n(r jeoes <n 
° o 1 


.«y guided rays; whereas, taking into account also 
sky rays [6]. this relation can be generalized 


2 Re « at wk ste 
n’(r }(eos*6 az sin?e sin?g ) ¢ ni. 


From the ray equations, one can obtein 





> r 
or e ee n*(r) ™ i?) * 12 4 
iw 35* mir) 0870, (—t90 etn) 1. 





so r{z) renges from 4 minimum velue R_ to 6 maxi- 
mum value R,° end is periodic with zr. ¢,° e,) 


as period [7]. R, and R, can be numerically obtei- 


ned letting equel to zero eq. (1). Then . one can 
get 2(rs ¢,° ®.) by the reletion 


R, 
1 or 
2) t2r.¢,0)-/ 2. 
2 0 0 0 4 
1 dz 


. f or 
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“8 ir} - (“tge sing 12-14. 
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For meridional rays, only R. exists, then one 
must let R,°0. For helicel reys, it is RoR, 209 
6lso in this case it is — 0. Finaluy. 
for the exiel rey. it is RRBs also in this ca 
se it is Zir_.@ .@ )*0. 6{z) end ¢(z) ere relae- 

oo 6 
ted to riz) by these equations 


n(r_) 





(3) qaeée(z) = cos9 |, 


n(r(z)) ° 


r te. eino. 
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therefore also 6(z) and @(z) ere periodic in z 
with 7(r)-@,°8,,) @s period. From the above equé 


‘fons, it is possible to see thet, for ref, or 
t°# . @ reaches its maximum or minimum value; 
@ 
⸗ 
in both these points. there is ¢* 2 . For helical 
, 
tays, it is uniformly 6({z)- —8 end ¢(z)-6° 3° 


for meridionel rays. it is ${z)-0." . 
Now we consider 6 well known class of index 
profiles. called a-profile, which hes the form 


n{r) = nf - 24( — ) « 


n? aa r? 





ere it is 4&* . By the zero-order WKS 


2n2 
0 


approximation [8], it has been demonstrated that 
the optimum velue of a in order to minimize the 
temporal window width of relative deleys results 


=)- — *1. 
ont 2-24. For this profile, considering a, 5 


end 4*0.02, by the numerical calculation of inte- 
grel (2). we have obteined the results which ere 
reported in fig. 1. 

















Fig. 1 - Curves of > 7(r.@, 6) vs. @, withr 
2 ° o 0 0 0 


es peremeter,. for an a-profile, with 


arta ef 1.5, 4*0.02. 
opt 0 


1 
These curves represents 3 


with r as peremeter. Only two different velues of 
. ere considered: ¢,°0 (dashed line), i.e. meri- 
& 
7 


Zir,@¢.0)vse.@, 
o oꝰ ꝰo 0 


dioneal rays: °° (continuous lines), i.e. skew 
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reys having rk, or roe - Precisely. the points 


2 


in which it results Z . O correspond to helical 


rays, the points on the left nands to reys having 
riko the points on the right hands to rays having 


8, Remark thet for a tending to zero, helical 


réys énd meridional rays have 4 po‘nt of junction, 

in origin of the plane (@ . 4), which represent 
the exial rey. Repeating these calculations for dif 
ferent valves of a, in the neighbourhood of a 


tre 


opt 
(from 2-44 to 2), we heve obteined no appreciable 
differences. 

Helical rey redii r.. ond the uniform propagation 
angle a, of such reys are related by this equation 


26 r So. 2 
oos ¥ 1° “Sein op 67 OP) rr, 


which can be solved numerically [9]. in fig. 2, 

MH 

= Vee 

é 

case of a =a 
op 


6, is reported for some values of a. In the 


* this dependence is neariy linnar, 


it is superlinear, for o > 6 
opt 


fo e 
Fe<G., it 1 


Sublinear. 





(ee 














Fig. 2 - Normalized helicel ray redii vs. the u- 
niform helical propagation angle e. for 


different values of a, m21+5. 4*0.02. 








3 - RAY OPTICAL PATH DESCRIPTION 
7 
To ol tain > by eq. (2) is correct, but, in or- 
Ger to echieve 8 precision better than 1 pert per 


10® , like the one which is required’ for the com- 
putation of the relative delay. starting from the 
ray optical path over a half-period, it is neces- 
sary to follow @ more rapidly convergent numeri- 
cal technique. This mathod consists of rewriting 
eq. {i) as 





R 
- (pleg 





(4) = \- n?(r) 


8 sing )*. . 
n2(R poy) 4 dazed 


where it is 
n(r ) 





we, ° nik) Oe, 


0 (for meridinal rays) 
sing, ⸗ 


(for skew rays) 


In this way, sterting from the initial condition 
(R ¢. 0). it is possible to solve eq. (4), 


with Runzge-Kutte’s numerical method, and obtain 

r(z) from R, to 8 If the calculetion is arrested 
F 4 

when r(z) reaches the value R,- 3 can be evéluated 

with @ very high precision. 

Substituting in (3), also 6{z) and ¢(z) will be 
obtained. These calculation can be extended to the 
second half-period, taking into account that r(z) 
and 6(z) result peir z-functions, whereas ¢(z) - > 


results 4 dispair z-function. In fig. 3 some plots 
of <. 8. @ vS. the normalized axial coordinate 


4 eZ J — — are rerertcd, for afta . Meri- 


dionul reys are conoteo fy cent inusus lines. For 





—* 
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Fig. 3 - Plots of = 6. @. vs. the normalized a- 
xiel coordinate z « 2/2(R + 6.) where 
it is ¢,70 (meridional reys) or “5 ° 


for ata 


» mn *1.5, 490.02. 
opt o 


simplicity. only two families of skew rays ere ta- 
ken into eccount: the first has R 20.3440 (dashed 


lines), the second R,°0.6908 (dotted lines). Each 

ray is characterized by o particuler value of *8* 

For meridional rays. J hes been selected starting 
from 0°. For skew rays, J has been selected star 
ting from the helical ray value 6. that is _— 

and J 6° respectively. 


From fig. 3.a and fig. 3.b an interesting property 
can be drown. For rays heaving the same Re R, in- 


creases with increasing *8* but the angle 6, which 
correspond to R,» remains the angle of the helicel 
ray having rok, This proposition is true also if 


one exchanges suffices 1 and 2. In other words, for 
whichever initiel radius ‘o" or 8° after «4 

, respecti- 
vely coincides with the helicel ray radius Ty which 


half-period <, the final redius R, or R 


corresponds <o the initial angle —8 


This property is schematically shown in fig. 4 














the protection cf the optical gputh of e generic 


shew rey on the trensverse plane 


is limited with- 


in two circuler trajectories, which belong to the 


helical reys heving redii Ta 


ting —* and 8.2 


“R.. * 
ond c. Let 


the propagetion engles of these 


helicel rays. the skew rey will have 6. in R. and 


G4 in R., 


H2 1 


- The enalytical demonstration of this 


essumption, ich is valid only for — will 


be done in appendix 1, sterting from the rey eque- 


tions. 





Fig. 4 - Projs<tion of a skew ray optical path on 


the trensve:=* olene. 
From fig. 3.c, it is possible 
for rays having the same Rye the 


greater amplitude of oscillation 
for reys heaving the seme e.° the 


smaller amplitude of oscillation 


to notice tiat, 
greeter o. the 


in @(2). whereas, 
greate: Rr, the 


in (2). In par- 


ticuler, for meridional rays (R,*0), it is 9(Z)*0 


(2<0.5). = (20.5). 


4. RELA'IVE DELAY CALCULATION 


Knowing r(z) and @(z), it is possible to obtein 
the group velocil'y of each ray (R,. ¢.° 6.) v6. zZ, 


which results 


viz) = cos (2) 


—_— 
n(r(Z)) 
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where < is the free spnce velocity of light. These 
velot tics are reportes in fig. 5S. for all the reys 
whose optical peths wie shown in fig. 3. In this 
computation the effects of muterieal dispersion are 
neglected. lhe possi) (ity of taking into account 
such effects has bee proved. om! will be argument 
of ↄ further work [9]. 

Meridiona] reys heve the lergest devietion from 
the fundamental velocity ‘ (which belongs to the 

° 

axial rey); then, the greeter R, the smaller de- 


vietions. Helical reys hewe uniform velocities, 
which are different from > but. in the figure, 6)! 
these velocities result nearly superimposed. 





%, 








4. 
* 











Fig. S ~ ®iot of the group velocities of the rays. 
weiose optice] psths were pictured in fig.3 
vs. 2; meridional rays in continuous lines 
skew rays having £140.344 and 81,0.69 re- 


spectively in dashed lines and in dotted 
lines. 


Meridionel] end skew reys trevel fester near R, 


then neer R.: Both fig. 3 and fig. 5. in which 


ore hes been considered. report results which 


eppeer nearly indipendent of a, in the neighbour- 


hood of —8 On the contrery, such little diffe- 


rences ceuse 6 very critical dependence on of the 
relative delays, either as sbsolute tempore) width, 
or (end especislly) as differential dGeley distribu- 
tion. 

















The relative Gelay of the generic rey (R+0,-8,) 


cer be calculated with respect to the celay of the 
exiel roy. vy the relation 


4 J 
° 
wo Oye 97 © uf vig) oc 





where L is the fiber length. It can be useful to 


c 
define « normalized raletive Jelay es * 1 Tt. 
⸗ a, 
Fig. 6 shows the reletive deleys es @ function 
of 8. for a*2-74; meridional reys are in dashed i* 


ne. helical reys in dotted line. Each family of skew 
rays is represented by e continuous line, starting 
from the helical rey point, towards the higher va- 
lues of 8.- All the guided mys present normalized 


2 
delays which renge from - * to 0. On the contrery. 
leaky rays. whose @reoa is denoted |. heve F ran~ 


ping from 0 to about a?. The relative Gelay of 4 ge- 
reric skew ray lies between the meridional rey line 


fr < » wh 4 Genote re- 
except for °°, —E ere °. and 40 


eanectively the helical ray having t equal to the 
minimum value and equal to zero. Approximately in 








of] 
. es ee — Se 
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[p=] P } 
ee ' 34 
‘* : ‘ jf 
e008 

















Fig. 6 - Relative delays vs. 6. with ®. 48 peramne- 
ters a*2-2A,. 0° 1-5-8490 .025 meridionel 


ray Gelays in dashed line, helical ray 
Gelays in dotted line. 








f , 
the middle of the interval le * 0) inericionel 


rey line and helical ray line intersect. 
Fig. 7 shows the relative delays foi a*2 end 


a*2-44. In the first case, 411 the guised reys he- 
2 


va 8 renging, from 0 to * +» On the contrary. lee- 


2 
ky reys h\eve 8 ranging from * t well abeowe a*. 


The relative delay of 4 generic sew ray is always 
limited by the helical ray line (/werds the higher 
velues of t. end by the meric! vl rey line towerds 
the lower values of x. The ca. @ of a*2-44 is the 
most interesting one, in fect 1 ranges ‘rom - 2 


to 0, for both guided end leaky reys. The relative 
Gelay of @ generic skew rey lies between the meri- 
dione] rey line end the helical ray line, except 
for a little @,-intervel eround the point of in- 
tersection of these two curves. Remark that, in 
this case, the leaky rey eres . is limited by the 
meridional rey line towerds the lower velves of 6, 


2 
and not by the horizontal line 6a * 


All these results ere in perfect agreement with 
[10], in which only helical ray reletive deleyswere 
computed, and with [11], in which it hes been de- 
monstreted thet the optimum value of a, taking into 
account also leaky roys, is 2-44. Finally they ere 
in good agreement also with [8]. in which only gui- 
Ged modes were considered. In the latter refe-ence, 
@ normalized propagation constant 6 was teken into 
account for denoting each guided mode; this pere- 
meter is defined as 


4 2 
6+ it. * 
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Fig. 7 ~ Relative <eleys vs. ®, with R, es pereme- 


<ers @*2 and a*2-44, 3. 4*0.02; me- 


ridionel ray Gelays in dashed line, heli- 
cal ray delays in dotted live. 


where 6 is the propagetion constent and Kon is 


the on-axis wevenumber: so it is O«éc4. The model 
veriable 6 and the rey-opticel veriables r.6 are 
related by this equetion 


6 - & nirieose 


from which, considering eleo the first of (3). two 
main conclusions can be drown. 

Firstly. many modes (or reys) heve the seme 6, 
and hence the seme 6; secondly. « relation between 
4 and 6, is necessery Cependent on f.. Then, to 


compere the reletive delay distribution t (6) ob- 


tained in [6] es averaged value over el] the modes 
heaving the same 6. with the reletive deley distri- 
but ion t (R08) of fis. 6 and fig. 7 is not so 


clear. Nevertheless, for meridional reys, this sim- 
ple reisation holds 


where suffi« M <‘s the abbrevistion of “meridional*. 
In [6]. it nas been odteined 


eo 2. *e 7 
t (6) 3 46) (a*2-24) 


(a*2-44) 


Compering these curves with the meridions) rey li- 
nes of fig. 6 and fig. 7. in which one scale of 6, 


is reported. it is possible to obtain « perfect 
agreement. 


A new technique in order to describe rey optice! 
petha — to celculate rey reletive deleys hes been 
presented. Comperisors with other techniques heve 
shown thet our method presents « time-resolution 
which is sufficient for giving « complete descri- 
ption of the model dispersion of fiberswith neer- 
“optimum index profile. Moreover, t’4 intrinsic 
tern of veriebles r. ¢. 6. which he) been used. gi- 
ves evident end intuitive optical «th pictures, 
which eppeer very interesting eapeci.\lly for dif- 
ferentie] mode measurements. In perti.vler, the 
effect of leaky reys hes been pointed owt. 
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From the rey equations, the expression r tee sine 


results « constant of the motion [7]; teking into 
eccount the only points in which it results 5" 


thet is Paes this reletion holds 


2?’ 
R,tge, — 

By @ series expansion, errestec at the 2-th order 

terms, one cen rewrite this equetion in the form 


($) Re, * Re, 


» it can be written 


From fig. 2. for or8 


ah 
r* ‘, 
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were & ts @ conétent. Therefore. in the sitwetion 


shown b fig. 4. in which eR *r » R* . 
' £ i it is a” ne a it 

results 

” R 

1 Z 

~~. o rr“ 

6 

m id 


Thit relation coincides with eq. (5). for Pi 
af 2 fo] — — . 
a wl 
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[Bext] apsract A chart is presented which allows one to design with good precision 
microstrip circular structures having an infinite attenuation at a given frequency. 
Effective parameters of the structure are automaticelly taken into account im the 
4e@qn chart, A mew type of bias filter element for active microwave integrated 
circuits is then proposed as an application of the chart. Notable sis reduction 
amd much better Gesign accuracy in comparision with the usual bies filter elements 


are demonstrated. 
1. INTRODUCTION 


The increasing interest in the study of 
microwave planar networks, i.e, of structu- 
rea with two dimensions comparable with the 
vevelength, is substantially due to two rea- 
coms, Firetly because, at high frequencies, 
‘te usual microwave integrated components, 
‘“eignued according to the standard monodi- 
thenmronal theories, exhibit behaviors unpre- 
‘.ttyble on the basis of such theories, Se- 
‘omay, because of the peculiar filtering pro 
mertis and of the lower impeawe?r levels 
presented by planar structures, 

It should be noted, however, that in 
*pite of the great number of theoretical and 
‘op Timental enalyees of such structures, 
‘hee design is still based on empirical 
tect niques, since the classical synthesis 
t+ tyods cannot be applied to them and re- 
t*mt attempts seem to be aenuselul, because 
uf oe — assumption of edge magnetic 
wale fi], 

At present, notable attention is devoted 
be planer structures with circular shape 
“+ 4), which, besides being used in the 
Constrection of circulators, present some 





(») art of this work was presented at the 6th 


Colloquium on Microwave Communication. 


Hudapest 1978. 


for each resonant mode, a good quantitative 
characterization of the behaviours of circu- 
ler microstrips, in such & way a8 to over- 
come rather simply the drawbacks yet un- 
derlined of other types of analysis. 

Starting from such # method of ana- 
lysis, a chart is presented in this paper 
which allows the design of microstrip twe- 
port circular structures having an imfinite 
attenuation at a given frequeccy. Effective 
perameters of the structure are automati- 
cally takén into account in the design char’, 
so that previous calculations of them arene 
required, A typical wee of the chart fro s#«- 
signing a bias filter clement is then ilius'tre 
ted and high cesign accuracy is demonsire- 
ted. A comparison with the usual rectencu- 
lar bias filter element shows a notable size 
reduction together with a comparable stop- 
band width. 

















27. 2H CINCULAR PLANAR NETWORK 


ihe filtering properties of the circular 
structure have been yet illustrated and explai 
ned both theoretically and experimestally ~ 
[5-6]; let as briefly report the cssential re- 
sults of its analysis. Fig. 1 shows the geo- 
metry of a microstrip circular structure, A 
Circular eection of radius 5 is connected 


with two feeding lines of width w; Wis the 
angle between the two lines and @ is the 


CS 


Fig. | - Geometry of the circular planar 
network, 


angle subtended by w. The normalized impe- 
dance parameters of this structure may be 
written: 
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where 
A — ẽ⸗ ff F om 
mn m eff mn fan 


4 ,- |i forms G 


x 2 form f 0 is the Newmanc factor 


oe * effective radias of the strecture 


J * dynamic permittivity of the TM 
_ mode [7] —_ 


* off * effective width of the input and out- 
put lines 


Fin * n-th cro of the derivative of the 
Bessel function of order m 
4 . 
eff a Cee “eff 
9 * is the vacuum light velo- 


It has been demonstrated [5] that pla- 
har networks present two types of trans- 
mission teros: modal and interaction ceros. 
Modal transmission zeros are located at the 
resonant frequencies of TM modes such 
that _ 


(2) cos’ im 1 


while interaction transmission zeros are 
located between the resonant frequencies of 


two consecutive TM and TM 
m 0 


i* ™,*° 


modes such that 


(3) agn(cos mn ¥) * sgnicos m, ¥). 


Therefore, keeping fixed the radius of the 
circular structure, the location of a modal 
zero is fixed, while that of an interation 
zero can be controlled by varying the posi 
tion of the ports, i.e. the angle (fig. 2). 
The first resonant modes of the circular 
structure are, in the order, the TM 509° 


the TM, 160 and the TM,, 9: It is worth ob- 


serving that while the TM... mode, which 


resonates at rero frequency, never gives 
clace to a modal tranemission sero, the 
oher two modes generally give plece to mo- 
dal tranemission sero, except for those va- 
lues of ¥ for which (2) is not satisfied, On 
the other hand, according to condition (3). 
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Fig. 2 - Frequency behaviours of the scat- 
tering parameter sj of the stru- 


cture in fig. 1, for different values 
of the angle y. 


the firet interaction zero is located between 


500 * 9 992 16 for v<90", and between 


° ° 
©. 50 and @,,5 fr 90< vw <180. it is 


easily seen, therefore, that in order to loca- 
te a transmission zero at a given frequency, 
notable size reduction of the structure can be 
obtained by choosing ¥<90. This is the 
case we have considered in deriving the de- 
sign chart. 


3. THE DESIGN CHART 
The exact location ®. of the interaction 


zero between © 500 * 0 and ®, 10 can be cal- 


culated by imposing the condition 
259 (@) +0 
which Leeomes, because of (1) 


a 2 cos mY 





a) F = ; +0 
m ono . w 
(g" nn! ett . . 2 ta'®! 
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Fig. 3 - The dew gn chart 


Let us define 


® 16 > ATL TT, _ 


as the frequency of the interaction sero nor 
malized with respect to the resonant fre. ~ 
quency of the TY 16 mode; eqn. (4) may 


then be writte, 











sr ¥ Ann co mY 


Mem alts) = © Hang), 9) 





It is worth noting that the coefficients aA” 
mn 
depend on the radius of the circular stru- 


cture and on the width of the ports, and that 
the ratio © ran'%)/ #1 (@, 19) is a function of 
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r o/h. where h is the substrate's thickness. 


As a coasequence, it is apparent from (5) 
that, for uw given substrate i.e. for given 
f and h, a is a function not only of whut 
r ’ 


also of w and * It can be found, however, 
that, by suitably limiting the range of varia 
tion of 58 and w, @ may be considered with 


good approximation as a function only of the 
wigie yw. In fact, for any and w such that 


r > 4h and w <r/4 the maximum change 
o 


of a is +1,5% for w =30 and less + 0.5% 
for y 245°. This leads to a substantial simpli 
fication in the construction and, consequently, 
in the use of the design chart. In fact, the 
structure can be designed in two steps: firstly 
by choosing the radius * and, consequently, 
by deterrnining the resonant frequency ® 10° 
secondly, by finding the angle wy such as to 
locate the transmission zero at the frequency 
prescribed, Incidentally, let us observe that, 
in the limit wheny — 90°, the solution of 
eqn, (5) tends to unity, i.e, the interaction 
zero tends to the modal zero located at 


o* ® 10° 


The design chart is shown in fig. 3; it 
is divided into four quadrants. 
I quadrant, The adimensional parameter 


0 
z ' z 
$= §',¢ co Pere V1 Mio? 





is showr as a function of yw. The straight line 
‘+ V = 1,84] represents the modal trans 


mission zero, whose frequency location, .as 
previously noted, is independent of yw. The 
other curve is representative of the intera- 
ction zero and has been evaluated numeri- 
cally through eqn, (5). 

Il quadrant. The adimensional parameter 


E reff 
n° Daf h h V 81 110? 








is quoted on the abscissae, The relation bet- 
ween » and ¢ represents a straight line whose 
angular coefficient can be found graphically 
through the non-linear scale on the top, or 
by making use of the auxiliary linear scale 
on the left and of the hyperbola, The con- 


struction is indicated in the figure in an 

example which will be illustrated later on. 
Ill quadrant, These curves allow one 

to determine from ro/h the value of 7 for a 


given substrate's permittivity €.° 
IV quadrant. When 8 and yw has been 


found, these curves are to be used to deter 
mine the maximum allowable value of the _ 
lines'' widht w such that the input and output 
lines would not overlap (taking into account 
their effective widhts). 

Some experimental results are also in- 
dicated in fig. 3, showing the good accuracy 
of the method of design proposed. 

It is worth notusm that one degree of 
freedom is left to the designer sii.ce Uither 
r, or y may be fixed arbitrarily, while the 


other parameter is determined by the condi- 


tion on the transmission zero, 
Let us describe in an example the de- 
Sign procedure, Sunpose 4* >= 4GI', .. ° 


= 10 (alumina), h = 0.0635 cm, w *: 0,062 cm 
(50 Q lines); let us fix r70.5 cm, Sinee 


r/h = 7.9, the value 9 * 25 is found 
through the curve €. = 10 in the Lil qua- 


drant, The hyperbola in the II quadrant is 
then used to find, starting from {jh ⸗ 


= 0.25 GHz cm on the linear scale, the 
corresponding point on the upper scale; the 
straight line connecting this point to the 
origin is the locus correspondent to the pairs 
rm ¥ which give place to an interaction ze- 


ro at a = 4GHz,. The intersection of this li- 


ne with 9+ 25 gives ¢ = 1.33 and then, 
through the curve in the I quadreat, w * 60°. 
Finally, the curve in the IV quadrant for 
w/h = 1 indicates that the value wp * 60° 

is acceptable, The experimental behaviour 
of the scattering parameter Is, vs, the 


frequency is shown in fig. 4. 
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Experimental frequency behaviour 
of the scattering parameter 18,01 


rig. 4 - 


of a circular microstrip structure 
designed to have a pole of atte- 
nuation at - = 4 GHz, 


4. THE CIRCULAR STRUCTURE AS BIAS 
FILTER ELEMENT 


In an active microwave integrated cir- 
cuit, the biasing network is generally obta- 
ined through a filtering structure like that 
rapresented in fig. 5a [8]. The double-stub 
section is normally designed according to the 


— 
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— — 








* 


Geometries of the usual bias filter 
element (a) and of the proposed 
new one (b). 


Fig. 5 - 


monodimensional theory in such a way as to 
present a short circuit at the frequency * 


of the r.f. signal travelling on the main arm; 


through a quarter wave transformer, the 


short circuit is transferred as an open cir- 
cuit at the main arm. The stop-band of the 
shunt arm can be increased by lowering the 
characteristic impedance of the double stub, 
i.e, by increasing its width 1, 

It should be observed, however, that 
for the usual values of such an impedance, 
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Fig. 6 - Scattering parameter Is.) vs. the 


frequency of the double stub of the 
BFE in fig. 5a, for two valucs of 

the stub impedance, In both cases 

the lenght of the stub is A/4 at 

fi, * 4 GHz, according to the muno 


dimensional design criteria, 


the stub section is, in practice, a two-di 
mensional structure and its behaviour 
cannot be predicted adequately by the clas- 
sical line theory. According to the two-di 
mensional theory, in fact, the pole of atte- 
nuation presented by the double stub, re- 
garded as a rectangular structure, is due 
to the interaction between the two consecu- 


tive resonant modes TMoo and TMo, and its 


frequency location strongly depends on the 
ratio 1/b between the dimensions of the rce- 
Ctangle. As a co: sequence, the characteri- 
stic impedance of the stub cannot be lowe- 
red too much, otherwhise the pole of atte- 
nuation would be shifted significantly from 
8 (see fig. 6 curve a), or even would not 


take place (fig. 6 curve b). 
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secause of the above considerations, 
it follows that, following the usual monodi- 
mensional design techniques, the transnius- 
sion zero frequency Cannot be predicted accu 
rately. unless a narrow stop band is acce- 
ptable. On the contrary, the good precision 
obtained in dimensioning a circular structure 
having a transmission zero at a given fre- 
quency can be usefully employed in the de- 
sign on a new type of bias filetr element; 
the above mentioned degree of freedom in 
the design of a circular structure may be 
used to optimise the stop-band width of the 
bias filter element. 

Fig. 5b shows the geometry of a BFE 
employing a circular structure instead of a 
rectangular one, The main line is supposed 
to have a characteristic impedance of 502, 
while the shunt branch have a characteristic 
impedance of 100 Q, 

By using the chart of fig. 3, several 
circular structures have been designed in 
order to optimise the stop-band width per 
cent of the BFE, The 30dB stop-band width 
% is shown in fig. 7 as a function of the 
angle w, for there values of the transmis- 
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ou — | » «gid = 4 4 
20 20 40 $0 60 70 60 v 90 
hig. 7 - 30d18 - stop-band width per cent of 
circular structures as functions of 
the angle y, with ‘. aS a para- 
meter, 
ston sere frequency te i.e, Of the signal fre 
quency, As can be seen, the maximum band 


sidth “s is reached for p between 40° and 60°, 
In the design of a BEE, the angle yw between 
the ports of the circular structure may be 
then determined by the condition of having 
the maximum stop-band width, 

The maximu:n stop-band width % is 
shown in fig, 8 as a function of the signal 
frequency fi similarly to the case of the 
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rectangle, this quantity decreases with increa 


sing * this is due to the fact that, by in- zs 
creasing i the structure's dimension decrea 
ses and, Correspondently, its input impedance 
level increases. 
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Fig. 8 - Maximum 30dB - stop-band width 


as a function of fo 


In order to complete the design of the 
BFE, the legth 1' of the line connecting the 
main arm to the circular structure should be 
chosen in such a way that an open circuit is 
transferred to the main arm for f * {5 In 


fact, since at the transmission zero frequen- 
cy fy the input impedance of the circular stru 


cture is low but different from zero, such a 
length should be different from A/4, (This is 
true also for the rectangular structure, al- 
through a study in this sense has not yet been 
published), 

Because of the fringing fields, however, 
the main line width and the circle radius are 
increased of the quantities dw, and ar, re- 


spectivelly in such a way that the phisical 
length 1' must be obtained by the relation 


vs Aw + de + 
2 0 


l 


where 1'', normalized to the wavelength A , is 
plotted in fig. 9 as a function of 8 fur those 


values of y which give place to the widest 
stop-band widths, 

As an example, the VSWR on the main 
arm obtained with a circular BFE with . 24 


GHz and v 40° is shown in fig. 10. Finally 
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Normalized length 1"' of the new 
type of BFE (fig. 5b) as a function 
of the frequency {5° for three va- 


Fig. 9 - 


lues of the angle wy. 
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Fig. 10 - VSWR on the main arm of a cir- 
cular BFE with w= 40° and for 
8 = 4 Gilz, 


the geometry of this BFE is compared in 
fig. 11 with the corresponding rectangular 
one, which has been designed according to 


the two-dimensional theory to give a trans- 
mission zero at the same frequency t. 24 


GHz. The performance of the two filters are 
quite analogous (stop-band width 21% and 19% 
respectively) while the former allows a nota- 


2 
ble size reduction (1.80 cm instead of 2,95 


em’). 
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Fig. 11 - Geometries of a rectangular BFE 
and of a circular one having equi 
valent electrical behaviours; both 
have been designed according to 
the two-dimensional theory [5]. 
Dimensions are in cm, 


5. CONCLUSIONS 


Starting from the two-dimensional analy 
sis of a planar microwave network with cir- 
cular shape, a chart has been derived which 
allows one to design with good precision a 
structure having an infinite attenuation at a 
given frequency. The simplest application of 
such a chart is in the design of a new type 
of bias filter element for active microwave 
integrated circuits. The limits of the usual 
design procedures for these filter elements 
are illustrated and are shown to be overcome 
by using circular structures designed by 
means of the present chart. It is shown more 
over that notable size reduction is obtained 
with the new type of bias filter element. More 
in general, the circular structure may be the 
refore usefully employed in complex circuits, 
whose elements have to occupy an area as 
small as possible. 
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EFFECT OF A TYPE OF MODULATION PULSE SHAPING IN FSK TRANSMISSION SYSTEMS 
WITH LIMITER-DISCRIMINATOR DETECTION 
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Consulente Telettra S.p.A., Vimercate (Milano)7 


[Text] Abstract. In this work we study a binary or multilevel continuous-phase FSK system with lim- 
iter-discriminator detection, in which the modulation pulses are rectangular with length 1 admit- 
ting of any value not higher than the symbol time 7. As an application example of the general 
treatment, we report numerical results in the binary case which, according to expectations, 
show the considerable inlluence on system performance, other conditions being fixed. In par- 
ticular, in the best operation condition of the system for P, = 10~* that takes place for the 
non-critical value r = 0.6 T, the transmission channel frequency band and the transmitting 
power are both smaller than those corresponding to the analogous best condition with the 
constraint + = T. Besides, the transmitting power advantages which derive from the opportune choice 
of ¢ increase little by little with channel frequency band reduction. With reference to multilevel sys 
tems, we have then point out some research development lines which we think interesting. 





1. INTRODUCTION 


A very import? nt aim of present research work on digital 
tadio-relay systems is the specification of spectrum efficient an- 
gular modulation techniques which allow use of highly non-lin- 
ear transmitters too. One way to rech such an aim is to make 
reference to systems which guarantee instantaneous phase con- 
tinvity of the transinitted signal se as to prevent it from having 
step changes. These, in fact, due to band limitation of circuits 
after the modulator, determine undesired signal arnplitude modu- 
lation which can cause great quality loss if non-linear power am- 
plifiers are used. 

This problem, which has opened up a very interesting re- 
search field, also mduces carcful consideration of the continuous 
-phase FSK modulation system whose simplicity, flexibility and 
reliability characteristics are well-known, especially when the de- 
modulator is a limiter-discruninator. To the authors’ knowledge 
the study of such a system has not been probed sufficiently to 
evaluate its performance im the general conditions imposed by 
the problem mentioned above. This performance has, in fact, 
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been cv lvated with teference to byt tanaainssion channels 
and for different cisaodulstors, taku; into account thermal 
noise and, a. must, intersymbol mterference or multipath fad- 
ing (sce for instance, from [1] ve [134). 

M: seover, Comparative evaluations of system quality with 
different shaping of modulation pul.cs do not appear te be a- 
vailable. As pulse shaping influcnces the modelated signal spec- 
trum, [34], 115], we think at as ompertant to take it mio ac 
count in a general analysis in which one of the bonds is the effi 
cient use of the frequemy band. This is what is done here where 


we analyse a binary of multilevel continuous-phase FSK system 
with luniter-discriminator detection, in which the modulation 
pulses are rectangular with length + admitting of any value not 
higher than the symboi time T. Our analysis generalizes the pre- 
vious ones, in which the modulation pulses have been consi4- 
cred rectangular with length r equal to T. 

The choice made for modulation pulse shaping has been 
suggested not only because of its particular simplicity but prin- 
cipally duc to the promising results obtained by applying the 
typical methods of the statistical decision theory to the corre- 
sponding continuous-phase FSK waverforms and shown in the 
next section. Then, with reference to the FSK system just spe- 
cified, we have expressed the c:ror probability P, by taking into 
acount thermal noise and the intersymbol interference produced 
by the lincar transmission channel. The expression obtained 
shera's im the first place that in the multilevel systems the usual 
choice of equidistant modulation levels [4] is not the best since 
the upper. positive and/or negative levels are unfairly treated with 
respect to the lower ones. However, the best choice of levels, 
strictly correlated with the length + of the pulses, is the object 
of current research on our part. 

In this work, as an application example of the general 
teatinent, we report numerical results in the binary case by 
providing graphics which allow both evaluation of the  influ- 
cur on system performance, other conditions being fixed, and 
tive best design of the system for P, = 10 . These results, 
which we wall cypose later on, show the considerable T influ. 
cnr, m accordance with previously mentioned expectations. 
In particular, in the best Operation Condiuiun of the system 
fr 1 10 © which takes place for the non critical value 


r 0.6 1, the transmission channel frequency band and the 
tranvnitting power are both smaller than those corresponding 
ty the analogous best condition with the constramt foi= 7. 
The transmitting power advantages which deriwe from the op- 
peortune chonwe of then intrease little by littl with channel 
hequetm ¥ vatd reduction, 

Finally, #¢ us interesting to observe that, whatever the value 
of rf (also t+ 0), possible errors in the sampling instants do 
hot m practice inctease the penalty with respect co the r = T 


iat, 
2 CONSIDERATIONS ON THE CHOICE OF MODULATION 
PULSL SHAPING 


Let us denote with 5, (t) and s(t) two possible waveforms 
defined m the interval (0, T,) having the same caergy 
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. To T> 
(1) =| shit) dt = | 53 (0) de 
0 0 


and the same a priori protabiluy P. Let us then consider the 
signal v(t) equal to s,(t) ov s,(t) plus a white gaussian nowse 
N(t) having two-sided power spectrum Ny tm the frequency 
domain. It is known [16] that the minimum error probability 
P., in the identification of s,(t) or s.(t) im eft) is expressed by 


1 w 
(2) P., = ~ exe Y= (1 —A), 





where 
E 
(3) 73 = 
is the intrinsic signal to noise ratio and 


1 /e 
(4) A= 7 5,(t) £,(¢) de 


is the normalized cross-correlation of s,(t) and s,(t). Because 
the condition 


(5) -1<A<1 


holds, (2) shows that, for a fixed value of W, P_. is minimum 
when A = — 1, Le. when 


(6) s(t) = —s,(¢); 


correspondingly we have 


(7) P —— erfe VW. 


Let us now refer to the case in which 
(8) s,(t) = Vy cos [(w, -w,) t+ 9), OS tS T, 
(9) s,(t) = Vy cos [lw +w,) t+ 0], Ot CT, 


with w, sufficiently greater than w, and 1/7, to practically 
verify (1). In such a case (4) gives 


by neglecting, as is allowed, the contribution of the s,(t) s5(¢) 
product at 2u,.. For a fixed value of W, the minimum of the 

















ert probsbidany PL. expressed by (2) takes place when A is 
minitnum, Le. from (10), when 


(1) wT, = 0.72». 


Se, the minana of A and PL. are 


(12) A= - 0.22, 


(13) 


Po, ee VOeiw 

2 

A binary cootmuous phase FSK siensl, with bic rate B, ~ UT, 
and with rectangular modulation pulxs having length 1 = T, 
comasts of a waveform succession of the type defined by (8) 
and (9. Thus, for this case, relation Sip (13) gwes a very import- 
ant reference: of Ws fined, P__ is preater than the corresponding 
minnwom error probability, because of the correlation among the 
waveforms. 

We should like to know whether an improvement of this refer- 
ence ts possible with the same W, by shaping the modulation pul 
ses in the PSK signal. A case easy to examine is that in which these 
pulses are rectangular agsin, but have length 1 admitting of any 
value in the interval (0, T,). So, let 


Vy cos [(w, — w,)t +0], O<et<r, 


(4) 5 (= | 


V, 


oe cos [wet — wr + 9}, r<t<T,, 


V,, cos[(w, +w,)¢+0], O<t<r, 


(1s) s(t) = 


Vy, cosluot +w,7+0], r<e<T,, 


be the two corresponding FSK waveforms, with w, satisfying 
the same conditions imposed in (8) and (9). By — 
(14) and (15) in (4) we obtain 


™-? 
+ cos 2 wf, 
T, 


sin 2w,7 
wT, 
neglecting, as in (10), the contribution of the s, (¢) s,(¢) product 


at 2 «2,. From (16), which coincides with (10) for r = T,, we 
castly deduce that A is minisnum as regards 7, w, being fixed, 


Ae 








(16) 


when 
. 
(17) w= 7T 
provided 
’ 
(18) “a? 9, 
because r < 7, must be true. 
From (16) and (17) 
Tt oe 
ay Awe 1 4 — 
T, 20,7, 
follows and, correspondingly, from (2) 
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1 
(20) P= — eke y~ (2- 
2 2 


in fig. 1 we compare the graphics of \ as a function of w,T, 

deduced from (10) where r = T,, (continuous line), and from 
(19) where 7 and w, satisfy (17) and (18) (dashed line). In the 
first place, we observe that by varying ¢ with «w, we can obtain 


2 o,7, 
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Fig. 1 — Graphs of normalized cross<correlation ) as a function 
of wT, deduced from (10) where + = T, (continuous line) and 
from (19) where + = 9/2, and wT, > #/2 (dashed line). 


values of ) smaller than the minimum expressed by (12) cor- 
responding to the case r = T,; so, with the same W, we can 
obtain values of P_ smaller than that deducible by (13). Be. 
sides, we observe that A tends toward - 1 and (20) tends to- 
ward (7), which represents the optunum. when in (17) w, + © 
and, therefore, 1 + 0; this on the other hand, can be imme. 
diately obtained from expressions (14) and (15) of the signals 
£, (t) and £,(t) which, on the basis of (17), satisfy (6), when 
r~ 0. 

Therefore, the considered shaping of the modulation pulses 
of the binary FSK signal allows notable improvement in the 
reference P__ and consequently in the minimum error probabil. 
ity when we operate on the waveform succession. This result 
leads us to think that the improvement may regard not orly 
the minimum error probability, but also the one we have in the 
real transmission and demodulation conditions of an FSK sy. 
nal. The investigation of such a possibility is the object of the 
next sections, in which we consider binary and multilevel FSK 
systems with limiter-discriminato: detection in the presence of 
intersymbol interference and thermal nowe. 


3. THE TRANSMISSION SYSTEM 


The block diagram of the transmission system which we con- 
sider here is reported in fig. 2. At the input I an L-level digital 
signal f(t) with symbol rate B, = 1/T, obtained by codifying an 
ergodic binary signal with bit rate B, = 1/T,, is incoming, It is 
well known that the connection between B, and fi depends on 
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Fog 2 Bho ke duagram of the transmucmon cyctem. 


i digital signal wpput, Nit) = thermal noise; 

Mo | PSK mealulator, Fl = receiving filter 

b= transmitting filter; D- = frequency demodulator (limiter discriminator); 
C= transmission medium, R = = regenerator. 


the code used. As an example, if L = 2 (1 = 1,2,3,...) we 
have 8 ~ Bi /l; consequently the relation T = 1 T, between the 
symbol time, T, and bit time, T,, is valid. As another example, aoe 
for the AMI code we have B = B, and so T = T,. 

In any case the modulating signal {(¢) may be written 





(2) fit) -), b, gt —h7), 








where git) is the function (fig. 3a) — 
“> -y* q : 
2 2 
(22) go- an 
6, otherwue, 


and b, a random variable admitting of L different values sym- R 
metrically arranged around zero, whose maxunum value we un- 


pose as |, without generality loss. — 


Fo + 














Fig. 3 Raat ay ae ges Nadas 
the angular frequencies w,. + b, w, }. 
—— of the function g(t) (when g(t) # Owe trans- 
mit the angular frequency w,, ). 














At th: carpet oft the conventional froquency modulstor Af we 
have 
129 Ayn Uy, cos fot + ate), 
where 1, amd or, are the aimpliude and the angular frequency 


ot the canner: the phase deviation a(t) 0s obtained by 


(24) 


a(t) - w fit) =, J b, git — 47), 


in which w, represents the mavimum angular frequency deviation 
of st). 
If we metroduce the function (fig. 3b) 


T T 

1, for -— +r ere—, 

. 2 2 

(25) z.(0= 


0, otherwise, 


we casily recognize that the signal A) can be written 


(26) Vy bo ae-ar cos [(w, +b, o,)¢ +o, 1+ 


Yad, £.(t—hT) cos lot + 1, 


where w, and w, assure the continuity of the s(t) instantaneous 
hase empoced by (23) and (24) To this end the relationships 


(27) (ws, +b, @,) (-- he nr) +m" 
= Ws (-- tht) +% 4, 
(28) (w, + 6, w,) (-- +r hr] +*“~* 


T 
= Wy (-— trehr) +, 
must be verified and from them we obtain 


(29) i —— +b, wr, 


(30) es oe 


(29) and (30) allow evaluation of #, and y, in correspondence 
with any 6, succession, when #, = a (0) is arbitrarily fixed. 
Let us now denote with 


(31) tt) = Vit) cos [wy + H(t) 
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lated by a single pulse bw, git - 


the transmission channel response at s(t) 2nd with m(t) that of 
the receiving filter F_ a: nowe Nit) which, as before, we suppose 
gaussian and white with two.sided power spectrum Ny in the 
frequency domsin. The output w(t) of the channel is therefore 


(32) 
The demodulator 2 here considered 1s an wleal limater discrim 
inator; its output is therefore $1). 


The regenerator R is composed of a sampler, a voltage compa- 


(33) - 1,6, 8,..<. 


w(t) = 4,(¢) + mle) = Uf) cos [we + (2)). 


(, =% +m, M @..., 
where 7,, which depends on the tranvmmsion system character. 
istics, oust be sunally chosen. 

The voltage comparator has 1. = 1 decison thresholds whose 
levels depend on the system and wm partwular on the 1. values 
which che random varizole by sysumes. In any case the threshold 
levels must be suitably determined « pusternuri; we oll unpose 
on them the bond to delimit some mtervals (with finite of 
infinite width) inside which, among the L ¥(t_) posible re. 
sponses of the system, only one response to the sygnal s/t) mod. 
mT), « included. The pulse 
generator following the compar ates then grcs at rts omtput the 
L-level digital regenerated signal. the corresponding syinbul suc. 
cession b, (h=...,~ 1,0, 1, . .) dors not coincide whith the 
b, one at the input because of decision errors determned by 
nowe N(t) and transmission channel distortwn. So, to evaluate 
system performance determination of the error probability P, is 
necessary. Such determination is the object of the next section. 


4. EVALUATION OF THE ERROR PROBABILITY P, 


In this section we evaluate the error probability P of the 
considered FSK system by making use of the method of paper 
[8]. To this end, with reference tu expressions (31) and (32) 
of 5, (t) and w(t), let us begin by determining the functions 
Vit), ¥(0) and (1). These functions can be easily expressed 
by introducing the in-phase and quadrature components of 
s{t) and mit). Asa matter of fact, if we put 








(34) s(t) = wif) cos Gigt — git) sin wet, 
(35) ait! = x(t) cos wot — y(t) sin wt, 
we obtain 
(36) u(t) = [p(t) + x(f)] cos wot ~lgit) · Am sen wt 
and from (34) and (36) the relationships 
(37) Vit) = Vp? (4) + 97 (0), 
pit) qi) ~ git) pe) 
(38) Mi 


p(t) +@°(#) 














1 
* “a Tole) + x0F +1¢in + OF 


tpi) ¢ xl gts) + Hol -igle + ylOl [ple + xte 


wnmedstcly dewend. 
To determine the values of V(t) and ¥(1) it is therefore necessary 
to calculate p(t) and g(t). To this end, let us introduce the complex 


envelopes 
ig=vy on, 
(40) 
it) = Vig) ef = ple) + igi) 


of the signals s(t) and « (¢). 
Recalling (26), the complex envelope i,(t) may be expressed by 


(41) 4S yg late MAE ge TM 


= a i,(*). 


where the meaning of the function i,t) ts clear. 
Let us now put 


(42) k=) ito. 


where i,(t) © the transmission channel response at i ,(1), and let 
us introduce the Fourier transforms | .(w) ar 1 (w) of i (1) 
and i.(1). If we denote with # (uw) and He) « + eransfer func. 
tions of the transmetting (1°) and receiving (Fj equivalent low. 
pass filters, and with C(w) the analogous function concerning 


the transinmsion medrumn, we have 


(43) glo) = Mfr) Cl) Hw) 1, (w), 


where, by virwe of (41), 











(= — 

Vv t =] 
(44) tw) = — f — = 

. T w-bhw, Ff 

uw 


44 am (al 


a (7) 


gg h® + fT) wae) , 











Prom (44), (43) and (42) we can evaluate i (f) and, conse- 
quently, p(t) and g(t). the characteristics of transmussion system 
being fin ed. 

As cegards error probabiluy 1”, we remember that the proba- 
bility distr tbuticn functsm 
(45) ig) » Prob (><a) 


of the random variable @ (¢) a1 the instant ¢ is known [17] and is 
expressed by 


(46) in) = A, Q (Va, Va,) + 4,11 - Ve, Ve,)1, 


in which: 
— Q (a, 8) denotes the Marcum function defined by 


(47) Qa, 8) = [ ge + #7024 (0k) dt, 
teense first kind and of 
~ the constants A,,A,, 4, anda, ave expressed by (') 

“i. 


wf 


Pe cates 
(7 M0) 9? +0} -20,9 


‘(Hs 2, -) | 
w} and o” are deducible from 


—— — 
"Wes $ug — 2u,9 











where w,. 


i - 
(50) “= 55/ w| H, (wo) |? des, 


i - 
($1) wi «= ~4 | w? | H,(w) |? dew, 


($2) o” = 2N, 8, 


im which 


1 - 
53 Be — 2 
(53) 9 fimo dw 


a ay 


The physical meaning of w,, w) and o” is clear; in partic- 
ular, o” ts che varience of the noise nif) ot 


frequency demodulator. 


at the input of Zhe 





It is opportune to mntrodyce im (49) the intrinsx signalto- 
poise ratio W which, assuming, for convenience, the insertion 
loss of the transmission channel equal to 1 at the angular fre- 
quency «,, id defined by (7) 


Vie 


We : ; 
aN, B, 





(54) 


So, the factor V" (4)/407 im (49) can be written 


win 3 BLM 
(55) (— W. 
4* 862 BN Wy! 


Knowledge of the function P(n) allows evaluation of the 
system error probability P_. As a matter of fact, let us make 
reference, as in [8], to the move interesting case in which ccw- 
sideration s allowed of sequences consisting of a finite number 
of symbols, and let us denote with S,, the symbol error proba- 
bility in the decision instant ¢_, for the k* sequence among 
those possible If the threshold levels corresponding to the 
symbol that the k®™ considered sequence presents in the instant 
t_ are 4, and 9), we have 

5,4 = Prob (4, <_)) + Prob (@, >2;) = 
($6) 
= Pin) + 1 — Pin). 


where 4, is the value of @(¢_) relative to the considered se- 
quence, for which the function P(n,) and P(qz) must be eval- 
vated by (46). 

The hit error probability 2, can be deduced from S,, when 
the conde & fixed. As an example, if L = 2? and we use a Gray 
code, the relation 


(57) r, SS 


is valul with very good appros uation. 


(8) te (4) ed (4%) the upper sign refers to the upper ity on the 
heft The quontitier a, peed a, ore prottiw on th Va, and Ve, iw (46) de 
fined a: thew prritwve square ruots. 

(7) We elwerve that in the binary case (54) and (3) coincide 


The crv predalabry " s then gwen by 


(38) PS EW), 


in which E(P,) os the expe. ted value of M,. 

The formulae we have obtamed alow ov dase of the mw 
trinsic ngnal toners tate Wowhah © mccoway om onder to 
guarantee a fixed error prob sbilay wlio the trawummon system 
characteristics are kftlown. For thes evalascan we hawe carried 
out 4 program which us a gencralization of that used mm [8] m the 
case L = 2! Up to this time the obtained results refer vo the bin. 
ary FSK system and are reporced in the nexi section. 

As regards the multileve! systems, we observe that the error 
probabihty depends both on the difference between the useful 
signal levels and the corresponding threshold ones, and on the 
values of the latter, as relationships (58), (57), (56), (46), (48), 
and (49) show. In particular, if we maintain such a difference 
constant, as has been generally supposed [4]. the upper postwe 
and/or negative levels of the signal are unfawly treated with 
respect to the lower ones. This can be easily understood in the 
very important case of a symmetrical transmission channel (ie. 
with w, = 0), by observing that in (49) the difference ¥ -- 9 
is dwided by a quantity which increases when| 9| increases 
Thersfore if the random variable b, is uniformly distributed, 
it is opportune to make reference to a chowe of the signal and 
threshold levels different from the previous one. Morcover, on 
the basis of the results obtained in the binary case, we think 
that also in multilevel systems the maxumum modulating level 
w,, strictly correlated with the length 1 of the mod: lation pul- 
ses, can be increased with respect to the case f = T, the trans 
mission channel bandwidth being fixed. However, on this sub 
ject the research is still in course of development. 


5. RESULTS OBTAINED FOR THE BINARY FSK SYSTEM 


Asan example of the previously obtained formu. 
lae, we have evaluated the performance of the considered FSK 
system in the binary case, with equiprobable and independent 
symbols b, = ¢ 1. So, the voltage comparator has only one 
threshold level, 9, whose optimum is q = 0. 

In particular, we have made reference to a system which is 
lacking in trasmission filter F,, employs an ideal transmitting 
medium C and has m the receiver a perfectly phase-equalized 
Butterworth filter F, with 4 poles. 

In order to evaluate its performance we have calculated in 
different cases the intrinsic signal-to nowe ratio W which is nec- 
essary to have an error probability P, = 10. *. 

As to the following it u convement to introduce the mormal. 
ined equivalent noise bandwidth b of the receiving filter, F, , 
defined by 


A. 


(59) 


_ and the normalized frequency deviation », connected to w, 
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by the relation 






















































































(60) —M 
ies, wies) 

In fig. 4 we report the values of W, as a function of duty , 
cycle r/T of the modulation pulses, necessary to have the fixed oj — — 
errot probability ,10* with the system considered, for * 
b ~ 0.9 and different values of #,. ee 

The optunurs sampling mstant has resulted in all cases im the 4 
centre of the modulation pulse; therefore, without lacking im 4 
generality, assuming null the phase and group delays of the fu ) 
ter, FL, we may write 

+6 T = — b-08 

Ve (oe) a-08 = —— 

‘is 12 = 
9 02 0.4 o6 068 i 1 
\ Fig 5 — Incrinsic signalto mows ratio W as a function of 1/T for 
1 _— ) P= 107°. for different values of b and for the values of o, which 
\ pove the optomum performance Notwe that along the curves the 
values of v, vary from point to point ; sce tables i, M1, M1, IV). 
- \ S3E5] The optimum for b = 1.1 and r/T = 1 is obtaiwd for v, = 0.34 
— \ / — and requires W = 12.8 dB. 
+s 
a e707 70 
* 22 o4 0.6 os 1 
v. b BS 

Fig. 4 — Intrinsic signaltonoise ra W as a function of duty 

cycle 0/1 of modulation pulses, for P, = 10° *. for some values 

vy mormeliced frequency deviation v, and for recewing filter 

eyumalent noise banlundth b = 09. 

Ter 

(61) Q2- > : 





Frown the curves of fig 4 we would first pomt out that for eve- 
ry #, we have noriceable variutions in system performance by va- 
rying the r/T catia, 

In porticu sat, every ¥, constant curve presents a minimum which 
takes piace for a value of 1/T which decreases as #, increases. More- 

over, we would pomt out that the value of 1 corresponding to the 
mintmem of every curve ms neat to the value #/2u, deducible 
from (17), for which the reference error probability reaches its 
minimum P, at least for relatwely high #,. 

in fg. $ we present fer sume values of bandwideh 6 the oyster 
perforrance obtamable by determining for every value of 1/T 
the frequency deviation », which corresponds to optimum, Le. 
which munemies W for P= 10° *. Then, along the 6 constant cur 
ves #, varies as a function of /T. Tables |, Ml, Ui, IV give the 
optumum values, ¥», of v, for different ratios 1/T and for the same 
values of b considered in fig. 5. In these tables we report also the 


quantities ¥, 1/T which, as may be seen, vary little by varying r. 











telle I 


Velurs of the imtrinee ngnalto-nowe ratio W, the normalized frequency deviation »,, and ¥,, 1/T which correspond to 


the optimum system performance for different values of +, for P, = 10 * and for normalized ceceiver bandwidth b = 1. 








































































































































































































ft Ooi; O1 02 03 04 os 06 07 os o9 I 
w 13.1 129 128 12.7 126 126 125 125 i25 125 126 
* 224 227 1.16 0.79 061 051 044 040 0.37 0.35 o™M 
¥_t/t 022 023 0.23 024 024 026 0.26 028 0.30 0 32 o™M 
Table li — The same as in table | except for b = 0.9. 
wt 0001; O1 02 03 04 05 06 0.7 08 09 1 
a 128 12.7 126 125 125 124 124 124 125 126 128 
* 226 2” 1.18 081 063 0.53 0 46 042 0.39 0.37 0.35 
¥,_tit 0.23 0.23 024 024 0.25 027 0.28 0.29 031 033 0.35 
Table Ill ~ The same as in table | except for b = 0.8. 
fT Ooo1; O01 02 03 04 os 06 07 os 09 i 
w 13 12.9 i28 128 128 128 12.9 13 132 13.4 13.7 
* 249 255 1.3) 0.99 0.70 0.59 051 0 46 0.42 039 03% 
*_0/t 07 026 026 027 028 0 30 031 0.32 0.34 0.35 0% 
Table 11° ~ The same as in table | except for b = 0.7. 
1 om; O01 02 03 04 os 06 07 04 09 
— — —— — — — — — — — 
“ 144 142 42 142 42 143 144 146 14.9 15.2 186 
* 2769 74 iS) 1.03 0.79 0.65 0.56 0s 045 042 039 
— — — — 
*,,*/1 0279 on” 0 30 0.3) 0.32 033 omM 0.35 0% 0.38 on” 

















liom ccammation of fig. 5 and the above tables we may 
make the tollowing remarks, 

The optimum for 1/1 = 1 is obtained forb=1 and 
vy O.38, an ches situation W = 12.6dit is required to guarantee 
e-08°. 

The optanam, considering alse rasa variable, » obtained 
teh OF F/T  O.6, Vay 0.46 and requires WW — 12.48; wah 
respect to the optimum Gase for 1/7 = 1, b and W are both sinaller, 
v, 1s greater. This helps to justify the last opinion of the previous 


section. The optimum is not critical: indeed the obtainable perform- (62) 


ance remains about the same for every 1/T in the interval 
(0.5, 0.7). 

When r ~ 0 the obtainable system performance varies very 
little with respect to the optimum corresponding to the consid- 
cred value of b if v, -+ . The optimum for r + 0 takes place 
for b = 0.9, v,7/T +0.23 and requires W = 12.8 dB. 

- The possibility of changing r enables us to get perform- 
ance improvements which are the more sensible the narrower 
the receiver filter bandwidth; for 6b = 0.7 we may see from 
the corresponding curve that it is possible to get an improve- 
ment of about 1.5 dB in comparison with the r/T = 1 case. 

Moreover, we have studied the sensitivity of the considered 
system to possible errors in the decision instant. Indeed, the 
choice cf narrower modulation pulses may lead one to think 
that greater sample instant accuracy is necessary. 
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Fig. 6 — Intrinsic signal-to-noise ratio W as a function of the 
sampling instat error, At/T, for P= 10 © for b = 0.9 and for 
different values of 1/T and v,. 


The results reported in fig. 6 show that such an opinion is 
groundless: the penalty due to an error At with respect to value 
r, expressed by (61), is rather small and fully comparable, in 
the cases 1/T = 0.6 (dashed line) and 7 very near to zero (dot- 
ted line), with that of the r/T = 1 case (solid line). In fig. 6 
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the curves are drawn only for At > 0, since for At < 0 they 
repeat themselves with parity law because of the perfect phase 
equalization of the considered transmission channel. 

Finally, for the optimum case, b = 0.9,1/f = 0.6, v, = 0.46, 
we report the graphs in one decision interval of normalized 
instantancous ampliude 1(¢)/1,, (fig. 7a), and of the noumal 
zed instantancous frequency deviation v (fig. 7b) defined by 


v 
2B. 





Because these graphs depend on the particular transmitted 
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Fig. 7 — Graphs of normalized amplitude V/V, (fig. a) and in- 
stantaneous normalized frequency deviation v (fig. b) in the case 
b = 0.9, r/T = 0.6, v, = 0.46 for the following two transmitted 
— ret ead emery pein mo 
—i-1, 1, ~1,—1,...(dotdashed line), 1 being the symbol 
corresponding to the considered decision eam a dashed 
line represents the normalized tr. ismitted frequency v git). 


symbol patternb, (h = ..,—1,0,1,...), we have chosen 
two particularly significant patterns. Namely, these are . . . , 1, 
i, 1, 1, 1,... (continuous 'ine) and...,-1,-—1, 1, —1, 
— 1, ... (dot-dashed line), | being the symbol c 

to the considered decision interval aa an ene 
T/2), without lacking in generality; the dashed line of fig. 7b 
represents the normalized transmitted frequency ¥,¢(¢). 





6. CONCLUSIONS 


We have determined the error probability P. of a binary or 
mltsleve | contmuous-phase 'SK system with hmiter-discrim- 
mator detection, in which the modulation pulses are rectan- 
gular with length 7 admitting of any value not higher than the 
symbol time T. As an application example, we have reported 
numerical results in the binary case which, according to expec- 
tations justified at the beginning of this work, show the consid- 
erable ¢ influence on system performance, other conditions 
being fixed. In particular, in the best operation condition of the 
system for P = 10° * that takes place for the nor-critical value 
+ = 0.6 T, the transmission channel frequency band and the 
transmitting power are both smaller than those corresponding 
to the analogous best condition with the constraint r = T. 

ve transmitting power advantages which derive from the op- 
portune choice of 7 then increase little by litthe with channel 
frequency band reduction. Moreover, we have verified that, what- 
ever the value uf 7 (also 1 -+ 0), possible errors in the sampling in- 
stants do not in practice increase the penalty, with respect to the 
rT ~ T case. Finally, with reference to multilevel systems, we have 
pointed out some research development lines which we think 
interesting and intend to examine closely in the future. 
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MSK, OFFSET QPSK MODULATION APPLIED TO A TDMA SATELLITE CHANNEL 


Milan ALTA FREQUENZA in English Jul 79 pp 284-289E 
[Rrticle by R. Agusti, Etsitb, Politecnico di Barcellona/ 


[ext] 


Abstract . In this paper we compare for the TDMA system conmuni- 
cation channel, the performance of two low envelope modulation 
systems, namely:MSK and OFFSET QPSK. Specifically, we consider 
a model of several cascade channels each of them consisting of 

a 4 pole Butterworth tranemitter filter, followed by a bandpass 
non linearity (TWT) and a Butterworth receiver filter. We su- 
prose that interchannel interference and thermal noise affects 
only the down-link., An optimal filtering method is found for 
the Gaussian interfering channels that achieve a given Bit Error 


Rate do with wsiniauve 5 ratio at the receiver input. The 


results are obtained by cobputer simulation. These results show 
a degradation of about 348 from the additive white noise ideal 
channel, when the normalized distance (D = fd T) among tranepon- 
ders is 1.333. 


1. INTRODUCTION 


The future dewand for high speed TDMA 
satellite links, points out to the need 
for modulation techniques, that saintai- 
ning the good efficiency of the PSK (Phase 
Shift Keying) signaling scheme, present 
low envelope fluctuations after filtering 
[1,2]. So, to obtain power efficiency, 
the power euplifier on-board (Traveling - 
- Wave Tube, TWT) @ust operate near seture 
tion, and consequently, in @ non linear so 
de. Otherwise, due to the unavoidable fil 
tering to reduce interference from the ad- 
jacent channels, the input signal envelope 
to the TWT is no longer constant and the 
Signal bandwidth spreads, producing sore 
Interchannel Interference (ICI). 


Once a sodulation technique is adopted, 
an adequate filtering is required to opti- 
mize the satellite channel performance, 
However, in the past [3,4], very little 
effort has been directed towards that end, 
because of the analytical difficulties 
with nonlinearities, even without ICI. A 
more realistic approach to @ satellite cha 
nnel must take this into account, 
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In this work, we have selected MSK (Mi- 
nimum Shift Keying) and OFFSET QPSK (Qua- 
drature Phase Shift Keying) as low envelo- 
pe sodulations, to present a Butterworth 


suboptiaal filtering configuration that 
sees sininal — . degradation for a given 


(10° 4) Bit mete Rate (8.E.R.). 
In fig. 1, we have shown three channels 


of our system, the one under consideration 
and the two adjacents ones, 
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Fig.l - Satellite Communication Link 


MOD Modulator 

TF Transmitter Filter 
DF Demultiplexer Filter 
MF Multiplexer Filter 
RF Receiver Filter 

DEM Dewsodulator 


Our hypoteses are the following: 


a) There is thermal noise and ICI only in 
the down link, That is, we assume large 


S ratio for the up-link and « DEMUX filter 


on-board that filters out the others cha- 
nnels. 


b) We treat the ICI as a Gaussian noise, 
to be added to the thermal noise, This 
hypotesis is useful to simplify the compry 


tation task and gives us an upper bound 
on the B.E.R. 


The system model that collects these 
hypothesis is shown in fig. 2, where we 
consider a Butterworth transmitter filter 
(TF) and a Butterworth receiver filter 
(RF) containing all the up-link and the 
down-link filtering respectively (the TF 
has 4 poles). 
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faterferance segecis 


é =H. 


Interference segneis 











Fig. 2 - Model for a simuletion of « 
digital satellite comsunice- 
tions link 


The received signal is heterodyned by 
reference sinusoids which are essumed to 
have correct frequency and phase. A Compu 
ter Simulation is used to obtain the BER. 


MSK and OFFSET OPSK have already been 
investigated by other authors [1,2]. In 
this paper we present some results for 
the more realistic channel previously men 
tioned that includ: an on-board non linear 
esplifier and ICI. 


Section 2 of this paper is devoted to 
a brief presentation of the modulation 
schemes investigated. In section 3 a 
known method for the computation of Power 
Spectrum after a nonlinear device is out- 
lined. Section 4 is the main part of the 
p‘per since Computer Simuletion and per- 
f<rmance evaluation methods are derived 
there. The final section contains the 
numerical results obtained in a particu- 
lar example. 


2, MSK AND OFFSET QPSK MODULATIONS, 

Both modulations are derived from cla- 
ssical QPSK modulation introducing a T/2 
delay in the data stream of the quadratu- 
re channel, in order to avoid excessive 
AM modulation after filtering. 


An OFFSET QPSK modulated signal is gi- 
ven by: 


e(t)=Re ([,Liarecty(t-nn) . jb, rect, (t- 


-hT- 5) |e Wot yone(B(t) © Jet) 
Where: 
Wo is the carrier frequency 
rect,(t/) = U(t) - Uft-T) 


(t) is the step function 
a(t) is the complex envelope of e/t) 





A MSK modulateé signal can be obtained 
with sinusoidal shaping of the pulses of 
OFFSET QPSK sodulated signal as: 


hoe a 
e(t)i=Re i 1 a, rect,,(t-hT)Singt+jb, 


co-@ 


rect,({t-hT~ F)coope] ett -pele(tier™**) 


After filtering, the modulated signal is 
given by: 


h== . 
r(t)=Re {ts p (t-hT)+Jb,p (t-hT- 7 


e¥et) = Re (r(t) eter) 


Where 


piti= 4 rectp(t)* k(t) in the OFFSET QPSK 


case 


pi(tj= i rect,,(t)Sin xtth(t) in the MSK 


case 


With het) the impulse response of the lov- 
-pass equivalent TF and r(t)the complex 
envelope of r(t) 


Assume now the pulse p(t) of finite du 
ration NT and define: 


how 
r v (t-h7,o,)e9¥"*) = Real 


h--@ 


ry(t)=Re { 


(Fy (t) „vet 
Where 
Py(t) is the complex envelope of r,{t) 


rh 


v(t-hT,¢,/= 9 mot) a,p(t-nT) + 


. b, p(t-nT- 5) |. recep (t-a7) 


Given that there are N tine slots of 
duration T whose content influence the 
real part of F(t), and N+) ones, influen 
cing the imapinary part of r(t), then 
v(t, ) is one of 22Ntl possible waveforas 


and the sequence {¢,) can be shown to he 


[5] a genes State Markov Chain, that beco 
mer homegencous and regular when{a.) and 
{b,) are statistically independent statio 
maty Sequences of independent binary ran- 
dom variables taking one of two equally 
likely values. 


The complex envelope of the modulated 
signal after the TWT can be expressed as: 


4 
a.(tJ= £ qlt-aT, o,) 
Characterizing the TWT as in [6], then 


s=10 
(1) q(t, o,)= L 


4 db, J, fealv(t, 1), 


v(t,o,) 


loct, 
Where 


b. are complex coifficients 
a is a scalar factor 


Ix) ie the first orden Bessel Function 
The modulated signal after the TWT is: 


a(ti= Re (3, (t) ef¥**) 


3. POWER SPECTRA 


Under the Gaussian hypotesis, the inter 
ference power is the variance of the Gau- 
ssian distribution assumed for ICI. There 
fore, we want to find the Power Spectrus 
of the modulated signal, ay(t), after the 
TWT. 


To obtain the Power Spectrum we have 
adopted a matrix formulation on ths basis 
of a L-state homogeneous Markov Chain 
whown in [1]. 


We summarize this method briefly. 
Let 


h-> @ 
y(tu= = 


h--e 


q(t - Af, ¢,) 
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Where @, is a L-state homogeneous Mar- 
kov Chain, with a transition probability 
matrix 


P = (pg4) ind = Uyecegh 
whose entries are 

Pes ⸗ Prob. (¢,=J/6, _,=t) independent 
of A. 


Provided that the Markov Chain is re- 
gular, i.e. that the limiting satrix 


P"= tie Pp” 
moo 


exits, the Pover Spectrus of &,(%) may be 
written as [1] 


Glw)=6(w) +6 ,fw) 


G (w= : Real {Q* (w) aA (w) Q@tw)) - 
- t o* fw) a (rT 4 P™) Qfw) 


. 2 + = — 2 
G 4fwd= a 2 (wiaP Qlw) J 5 (w-m 7 
m=-@ 


Where 


+ 
( ) means conjugate transponse 
I is aL * L identity matrix 
Q(w) is a column L-vector whose i-th en- 


try is the Pourier Transfores of 
q(t, ,/ 


” is -a diagonal matrix with entries 


fajij = 








P, = Prob, { +, = i} 


A(w) ie a matrix series 


Afw) = r (Pp = p™y"Q~ smut 
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A problem arises in evaluating A(fw). 
Although a computational algorithse has 
been devised [7], this become impractical 
when W>2 since we should need to store in 


the computer sesory 27m matrices each 
of thes having 22N+l 27+l elements, i.¢., 


2?! words if wNe-3. 


We2 gives us a good aproxisstion, 
d,(t), to a(t) provided that Beq,T(Beq, 


is the bandwidth equivalent noise of the 
TF) is large enough (practically, Beq ,? = 
1.5 is sufficient). If Beq,7? < 1.5, the 
sentioned algorithms is no ibnger useful 
since ve aust take N>2, 


A solution consists in adding only 
the first terms of A(w) provided that the 
Sequence (P") converges rapidly enough as 
m increases. Then, N#3, gives us a rea- 
listic approximation, #,(t) to #(t) using 


only # reasonable amount of semory. This, 
can be further reduced observing the exis 
tance of a delay in thé appearance of sig- 
nificants values in p(t) when filters are 
used. So, if we assume 


Tr 


p(t) = 0 for o<t< g 


a 732 reduction in the memory occupancy 
is gained, i @., P has nov 228. 22" ele- 
sents, 


A transition matrix that leads to «4 
rapid convergence is buit in the following 
way: 


We define 


bY, e142(a',+1)44(b',_ y¢1)+8(a',_;* 





1+ 
¢,= 2 


41)416(b',_s¢1)+32 a", _,*1) 
2 





b™,_ +142 (a",+1)44(b",_y¢1)* 
* 





", +1) 
+86(a",_,+1)416(b", _ +1) +32 (0 h-2 — 





Pe3* Prob. (¢, ,,=51¢,=t)= 3 if 





Pog = 0. In other case 


a's, a”. and b's, db”; are two saeuples of 
the reandows variables a, and b; respective- 
ly. 

4, SYSTEM PERFORMANCE. 

4.1, Caleulation of interference power 


Since q(t, @,) can take only an even 
nuaber of equally waveforms vith syasuetri- 
cal properties: 

Gz (wv) =0 
Then 
Glw)=6 (vw) 


And the Caussian variance is given by: 


2_1 ~ ¥ dv _ 
oy J. , Gl fin. fd) lle) Fe = 
a GC, tw) Btw) de 
4 2s 


Where 


The 1/2 fector is introduced to ob- 
tain the real power from the invelope po- 
wer. 


* 
Hiw) ie the Fourier Tranafore of A(t) 


f is the spacing frequency separa- 
tion among transponders. 
A Simpson quadrature rule is used to 
Compute 2 » The bound fy, required for 


the integration limite is settled accor- 
ding to the error specified. 


Error = fy Gy (wo) Atw) * 


From the computation results of Gy (w) 
we have G,(w)<2T 





Then e 
d 
ee e -2n 
Error<2T Jr, 7 aly crit dy 
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2feT 
(2n-1)(f,/ fo 


< 





fe ie the 3 4B bandwidth of Aw) 
n is the number of poles of A(w) 
f, has been selected to ssintain an Error 


< 107? 

When Beq ,=1.5 both algorithas shown 
in section 3, lead to the sane o* as vasto 
be expected with that intermediate velue. 
4.2. Computer simulation 


We perform a computer sisulation to 
find the bit error rate. 


The signal source provides the soduls 
ting signal ft) in both OFFSET QPSK and 
MSK cases. 


The 
ding to 
senples 


filtering is carried out accor- 
@ state variable method taking 20 
in T seconds. 


The TWT is simulated by (1) 


The demodulation process is perforsed 
by a perfect coherent deaodulator in the 
sense that if @¢— is the mean value of 
received phase, then the demodulated sig- 
nel, raft) is: 


ra(t) = 5 iy (t) » hee) eit 


Sampling the real and imaginary pert 
of raft) at the points t + AT, we obtein 


respectively the {a,) and {b,) sequences 
of detected data, 

In presence of Gaussian noise of sero 
mean and 6 variance, the BFR per channel 


can be expressed if the eye ciagrama is 
Open as: 





P 1 la’, | 
bp =f {z erfo—— } 
/20 
|b’, | 
] h 
P = E {gerfe } 
bq 7 e 
Where 
2 2 
— —50 I + NyBeq., 
Beq ie the bandwidth equivalent noise 
of RF 


Ne is the one-sided noine spectral den- 
sity of thermal noise. 








Thea the symbol rate is 


Pole) ia ceadily simulated when a aaxi- 
mal length PN sequence is assumed as modu 
lating data, i.e, if there are M slots 
producing intersyabol interference, Choo 
sing a 2 =- 1 PN sequence, all the possi- 
ble combinations of assumed interfering 
data are present, then 


nz! — 


erfc 


























p 24 1 r s 
bp 7 ‘eT 24 2 
M-1 , 
1 1 “yy t la’, 
*- 3 · * erfe 
z n=l eq,-T 2 
1 ° 
P = J — 7 erfc pa "| = 
bq z° a n=1 
1 1 nea , f et 
2 — aa erfc — 
7 n=1 = 2 
n 
' ; 
Where 
gE 
n= — » with Fy the bit energy prior 
Ne the TWT 
Pole) is computed for the 20 possible 


sampling points, and we select te as the 
best in the sense of sininal Pie). 


4.3, Optimization procedure, 


In order to perform the filt cing 
optimization method, we start «suming 
for the BeqeT of TF a value of 0.9 and 
incrementing it in steps of 0.1. For 
each chosen vajue we compute the inaterfe- 
rence power, 0 r for all the possible FR 
(n*poles, 2,3,4-55,5 and Beq 7 = 06.8, 9.9, 
Tyee) 


P fe) ia computed only for n = 12,13 
and 19dB, A graphic linear interpolation 
carrier out on double logarithsic paper 
(Well behaved for 1448 > 5 > 1048) allows 
us to determine the RF that introduces 
the least degradation in 4 for P. fej = 
= 2,107", Thies method is suggested by the 
asyaptotic formula. 
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2 
erfe (zx) ~ e* 


from which follows that 


log ( log 





} = 
erfc(z) 
= 20 log z + 10 log (log e) 


To get a best approximation, the ob- 
tained mn is recycled in the siaulation run 


until we reach Pi fe) — 2.10°° 
5. RESULTS AND CONCLUSIONS 


We have adopted a f, = 80 MH. and 
@ transmission rate of [20 Mbit/s. 


The values of n required to obtain 
the 107 ‘oer are represented for the MSK 
case in figure 3, in which the three best 
RF (specifyed by a number, see table i 
where we have listed the number assigned 
to each Butterworth filter) are considered 
for each Beq,T of the TF. 








Table 1 
a “poles Beq..T 
1,8,9 6 Be: 
2,3,4,5,6,7 4 i 
10, 11,12 6 1 
13,14,15 3 1 
16,17,18 3 0.9 





























"n 4 4 = A A. 4 i> 
’ u 12 1 4 1S Seq, 7 
Fig. 3 = versus Beq. T in a MSK systen 


with ¢,T = i*3333 





The siniaunm n is found to be 11.7 dB 
to which it corresponds a TF with a Beq,7 


~ 1.3 and ar with 4 poles and « Beq .f = 1 


In the same way figure 4 and table 2 
shows the system perforeance with an 
OFFSET OFPSK modulation. The best 4 is 
now 11.6 48, and to this value corresponds 
a RF with 4 poles and Beq .f * 1, as befo- 
re. 

Hovever, the optimum is not as well 
behaved as in the MSK case, since a flat- 
neas in the curve appears, due to the pr 
sence of less interference power (table 3) 
with OFFSET OPSK sodulation, when « tigh 
channel frequency spacing is adopted. Any 
way, @ small value of BDeq,? gives us « 
more realistic hypothesis of No ICI pri- 
or the TWT. So, feq,7T7!.3 can also now 
be adopted with nz11,7d8. 








Table 2 

n*po 

‘les Beq,T 
1,2 6 1,1 
3,4,5,6,7,8,9,10,11,12 4 i 
13, 4 1.1 
14, “ 1 
15,16,17,186,19,20,21,22,23,24 6 1.1 
25,26,27,28,29,30,31,32 
33,34,35,36 3 i 








— 
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. > @ * on 
4 bea, Tt 
ge ye — — — — — 
lee ee 2 ee ee ee es J 





Fig.4 - ) versus Beq tT in a Offset QPSK 


eyetees with £,.T = 1,333 


a4 











Table 3 
—2 
Beq,T a MSE 0° OFFSET OPsK 
" 0.015 0.013 
1.1 0.018 0.014 
1.2 0,021 0.016 
1.3 0.025 0.018 
1.4 0.030 0.020 
1.5 0.035 0.021 
1.6 0.041 0.02% 
1.7? 0.047 0.024 
10 0.089 0.050 
Beq tT - 1 


RF has 4 poles 





Therefore, both schemes of modulation 
present an identical optiwal system per- 
formance, i.e., the beat 4 remains sbout 
34B bellow the ideal case. 


A comperative analysis between this 
two sodulations scheses, according to the 
chennel frequency spacing used (fd), 
would be interesting. Wowever, [ros 
these results we can assume that if 


fal < ay = 1.338, 


the OFFSET QPSK scheme has better perfor~- 
mance then the MSK one, and converse: ly, 
when f,.7 > 1,433, because of the idcenti-+ 
cal cyften performance at i = 43,4387 in 


the optimal case, and the higher spectral 
density power in the neighborhood of car- 
rier frequency of the OFFSET QPSK signel. 
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